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UNIFORM EXPONENTIAL STABILITY AND CONTROL 
CONVERGENCE OF SEMI-DISCRETE SCHEME FOR A 
TIMOSHENKO BEAM * 


FU ZHENGİ, ZHEN JIA?, AND BAO-ZHU GUO 


Abstract. This paper considers numerical approximations of a Timoshenko beam under bound- 
ary control. The continuous system under boundary feedback is known to be exponentially stable. 
Firstly, the continuous system is transformed into an equivalent first-order port-Hamiltonian formu- 
lation. A basically order reduction finite difference scheme is applied to derive a family of semi- 
discretized systems. Secondly, a completely new method which is based on a mixed discrete observ- 
ability inequality involving final state observability and exact observability is developed to prove the 
uniform exponential stability of the discrete systems. More interestingly, the proof for the uniformly 
exponential stability of discrete systems is almost parallel to that of the continuous counterpart. 
Thirdly, the solutions of the semi-discretized systems are shown to be strongly convergent to the 
solution of the original system through Trotter-Kato theorem. Finally, both exact controllability 
of continuous system and the discrete systems are proved in light of Russell’s “controllability via 
stability” principle and the explicit controls are derived. Moreover, the discrete controls are shown 
in first time to be convergent to the continuous control by proposed approach. 
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1. Introduction. A big concern in control of systems described by the partial 
differential equations (PDEs) is the finite-dimensional approximation. This is because 
most of the feedback controls like observer-based feedback controls for PDEs are 
usually infinite-dimensional. In order to apply the designed PDE feedback control, 
the discretization is indispensable in most of the situations. For this purpose, a 
first natural step is to discretize the spatial variable of PDEs and at the same time 
to keep the time continuous. This process is called the semi-discretization. The 
benefits of such semi-discretization is obvious. On the one hand, the semi-discretized 
models are finite-dimensional which can be handled by most of the control engineers, 
and on the other hand, if the semi-discretized modes preserver as many as physical 
properties of the original PDEs, the discrete modes themselves can be regarded as 
physical modes instead of infinite-dimensional PDEs. Unfortunately, this is not always 
possible because the spurious high frequency modes might be produced during the 
discretization, which destroy such a preservation. The papers [2] and [I4] first pointed 
out independently that the classical finite-difference or finite element discretizations 
for a 1-d wave equation can not preserve the uniform exponential stability. This has 
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been confirmed by many researchers in latter two decades, which includes not only 
the uniform exponential stability but also the uniform controllability and observability 
[31] [38]. This seriously restricts the applications of the PDE feedback controls. 

Several remedies have been introduced to overcome this difficulty like those Ty- 
chonoff regularization [I4], mixed finite elements [3] [4], filtering of high frequencies 
[3] (15), two-grid algorithms [27], non-uniform numerical meshes [8], and vanishing 
viscosity approach [25] [26] BI]. Among many of these aforementioned methods, the 
vanishing viscosity approach is a popular one, which introduces a vanishing damping 
term with step size in discrete models. However, these approaches either lose the 
popularity of the finite difference scheme or add artificially numerical viscosity which 
is different for different systems. In addition, the artificial numerical viscosity method 
is most often confined to special boundary conditions because the proofs for unifor- 
m exponential stability and controllability rely heavily on the analytic forms of the 
eigen-pairs. Very recently, a natural finite-difference scheme based on order reduction 
method was introduced in 22] for wave equation with boundary damping. It has 
demonstrated that the scheme not only preserves the uniform exponential stability 
but also the proof of the uniform exponential stability is completely parallel to the 
continuous counterpart, which simplifies significantly the mathematical analysis. The 
approach is somehow universal, which has been applied successfully to uniform ex- 
ponential stability for wave equation with local viscosity damping ({13}[37]), uniform 
exponential stability of Schrödinger equation ({12}[20]), uniform exponential stability 
and observability of Euler-Bernoulli equation ({24]), and even observer-based uniform 
exponential stability of wave equation which is actually a coupled PDEs ([30]). The 
left problem however, is the convergence of the discrete controls, which has been out- 
standing for quite while. The uniform controllability of the semi-discretization for 
beam equation has been studied in [6] [7] by numerical viscosity approach. 

Motivated from these works aforementioned, in this paper, we study a classical 
Timoshenko beam which was extensively investigated in the past decades due to its 
widely applications in engineering, both from PDEs and numerical approximation 
point of view. The Timoshenko beam model is a variant of the Euler-Bernoulli beam 
model, which takes the effect of shear and rotation into account. In this paper, we 
consider a Timoshenko beam model described by the following equations: 


pwr (x,t) — K(w,(2, t) — (x, t))s =0, xz € (0,1), t>0, 

I, bu (2, t) — EI¢z,(2,t) — K(wz(x,t) = (2, t)) E 0, LE (0, 1), t> 0, 
(1.1) w(0,t) = ¢(0,t)=0, t> 0, 

we (1, t) = (1, t) = u(t), t>0, 

Px(1,t) = u2(t), t > 0, 


where w(x, t) is the transverse displacement and ¢(a, t) the rotation angle of a filament 
of the beam, both at time t and location x. The constant coefficients p, Ip, EJ and 
K are the mass per unit length, the rotary moment of inertia of a cross section, 
the product of Youngs modulus of elasticity and the moment of inertia of a cross 
section, and the shear modulus, respectively. ui(-) and ug(-) are boundary controls 
that monitor w(x, t) and (x,t) at x = 1 and transform them into the lateral force 
and moment applied at x = 1, respectively. For notational simplicity, hereafter we 
omit all obvious space and time domains in the rest of the paper. 

Stabilization for Timoshenko beam has been studied by many researchers through 
designing boundary or interior damping in the past several decades. Kim and Renardy 
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proposed the following boundary feedback controls in paper [I8]: 
(1.2) u(t) = —a pu;z(1,t), u2(t) = —ael,¢:(1, t). 


They proved that the energy of system (L.I) under boundary controls (1.2) decays 
exponentially by the Co-semigroup theory and the Lyapunov functional method. In 
[34], Xu and Feng presented detailed spectral analysis and Riesz basis property of 
the generalized eigenfunctions of the beam (L.i)-(L-2). Gorrec, Jacob, and Zwart et 
al. obtained the exponential stability of (LI)-(L2) by the stability theory of port- 
Hamiltonian systems on infinite-dimensional space in [33]. Raposo et al. added 
interior frictional dissipative terms w(x, t) and ¢;(zx,t) into (LI) to get 


pwr(x, t) — K(w,(ax, t) — (x, t))s + wlz, t) = 0, 

Ip bit (2, t) ~~ Elézz(a, t) ~~ K(wz(2,t) ~ olx, t)) = hilz, t) = 0, 
w(0,t) = ¢(0,t) = 0, 

wz (1, t) a (1, t) = bx(1,t) = 0. 


and obtained in [28] the exponential stability of system (1.3) by means of the frequen- 
cy domain approach of Co-semigroups. Yan et al. studied in [35] the stabilization 
problem of Timoshenko beam (I.I) under linear dissipative boundary feedback con- 
trols and presented various necessary and sufficient conditions for the system to be 
asymptotically stable, where the equivalence between the exponential stability and 
the asymptotic stability for the closed-loop system was also given by the frequency 
domain analysis. 

Recently, Rivera and Naso proposed boundary dissipation only on one side of the 
bending moment, i.e., EId,(0,¢) = —a2¢;(0,t) with other three Dirichlet boundary 
conditions to (LI). They proved the exponential stability under the condition that the 
wave speeds of the system are equal to each other [29]. Almeida Júnior, Ramos and 
Freitas found some new facts related to the classical Timoshenko system. Precisely, 
they proved the damped shear beam model, which corresponds to a part of the classical 
Timoshenko beam model, possesses an energy exponential decay for any coefficients 
of the system under the Dirichlet-Neumann boundary conditions [I]. 

In this paper, we are focusing on system (1.1)-(@.2) and apply the Russell’s prin- 
ciple to construct controls which steer any initial state to rest both in both continuous 
and discrete cases. However, it is well known that either the exponential stability or 
uniform exponential stability of some systems are needed when Russell’s principle is 
applied. We hence start from the uniform exponential stability of discrete systems 
since the exponential stability of continuous system has been well established. The 
most difficult part of this paper is to test uniform exponential stability of discretized 
Timoshenko beam, for which we have encountered difficulty in applying existing ap- 
proaches. The main reason is that it is difficult to find a suitable Lyapunov function 
or a frequency domain energy multiplier for this coupled system. In this paper, we 
propose a completely new proof for the uniform exponential stability, which is quite 
different from previous works [25]. The convergence analysis, i.e., the solutions of 
the semi-discretized systems strongly converge to the solution of the original system, 
is also developed. This is realized through Trotter-Kato theorem and do not utilize 
the methods of [20] [20] [3I] because the results there can not be applied directly to the 
discrete approximations of the controls. Based on the uniform exponential stability 
result and convergence analysis, we derive that the discrete controls are convergent 
to the continuous system control when the initial values of the continuous system and 
discrete systems satisfy some priori convergence conditions. 


(1.3) 
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The rest of this paper is organized as follows. In Section 2] we restate the ex- 
ponential stability of an equivalent order reduced continuous system and discretize 
its spatial variable by the method of [22] [20]. Uniform exponential stability is also 
presented in this section. In Section B] convergence of the numerical approximate 
solutions to the continuous one is analyzed by using the Trotter-Kato theorem. In 
Section [4] we obtain that the controls that steer the initial values to rest by Russell’s 
principle both for continuous and discrete systems. Moreover, we show that the dis- 
crete controls are strongly convergent to the continuous system counterpart, followed 
up by concluding remarks in Section 5] 


2. Preliminary results on systems (1.1)-(1.2). In this section, we discuss 
uniform exponential stability for Timoshenko beam (Li) under boundary feedbacks 
(7.2). Subsection 2.1] gives exponential stability for continuous system and Subsection 
[2.2] presents uniform exponential stability for semi-discrete models. 


2.1. Exponential stability of continuous system (I.1)-(1.2). To have an 
almost one-to-one correspondence from the exponential stability of the continuous 
PDE system to the uniform exponential stability of the discrete ODE counterparts 
discussed in next subsection, we introduce order reduction transforms: 


yi (a, t) = Wz (2, t) E olz, t), yo(z, t) = pwr(2,t), 
(2.1) 


y3(x, t) = bx (2, t), ya(x, t) = o9t(2,t), 


and obtain an equivalent system of (1.1)-(@.2) as follows: 


yi (a, t) = B2yo(x, t) E Baya(z, t), 
yo(a, t) = By, (2, t), 
(2 2) Yys(a, t) = Baya (a, t), 
l ya(x,t) = b3y3(x,t) + biyi (z, t), 
yi(1,t) = —-a14Y2 1,t), y2(0, t) = 0, 
ys(1, t) = —a2y4(1,t), y4 (0, t) =0, 


in which 6; = K, 62 = p~', 63 = EI and $4 = I>", and dot’ to stand in shorthand 
for the derivative in t and prime / for x. It is seem that in system all boundary 
derivatives with respect to x disappear, which is the biggest merit of the order reduc- 
tion method. Set Y(a,t) = (y1(z, t), y2(x, t), ya(x, t), ya(z, t))'. Then the first four 
equations of (2.2) can be written as a port-Hamiltonian system 


(2.3) Y (x,t) = P,[HY (2, t)]’ + PoHY (x,t), 
where 
0 0 0 -1 0 1 0 0 
, jooo 0 _{1 00 0 
H= diag(81, b2, 63, Ba), Po = 00 0 0 ; Pı = 000 1 
1 0 0 0 0 0 1 0 


We consider system (2.2) or (2:3) in the state space X := L?((0,1);C*) with the inner 
product 


(2.4) (YZ), = T Y(£)"HZ(z)dz, Y Y,Z € X. 
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The rest of this subsection has been employed in [33]. To formulate the boundary 
conditions of (2.2) or (2.3) in a better manner, we introduce the boundary effort and 
flow defined respectively as ea = vV 2-1H[Y (1) + Y(0)] and fa = V2-!P,H[Y(1) — 
Y(0)]. We therefore consider the operator 

AY (x) = PI(HY (x) + Po(HY(a2)), Y Y € D(A), 
(2.5) 


Baje T €X: HY € H1((0,1);C), Wg £] z o} , 


where H!((0, 1); C4) is the Sobolev space and 


a, b28} 0 0 BoBy* œ 0 0 
1 =f 0 0 0 0 1 0 (0) 

2.6 Wg = — os = 
ee) amo 0 as Babs? 0 0 fabs? ap 
0 0 —1 0 0 0 0 1 


In this way, system (2.2) is formulated as an evolution equation in X: 
(2.7) Y(-,t) = AY (t), Y(,0) €X. 


The following Lemma gives basic property of the matrix Wg. Because it involves 
simple calculation, we omit the proof. 
LEMMA 2.1. Let I, be the identity operator on Ct and £ = f 
4 
matrix Wg has full rank and satisfies Wg &YW} > 0. 
Next Lemma P.2] was brought from Lemma 7.2.1]. 


LEMMA 2.2. For any Y € D(A), there holds 


(2.8) Re(AY,Y)y = —a1b1b2ly2(1)|? — a2838alya(1)/?- 


Now we can show that (2.7) or (2.2) is exponentially stable. 
THEOREM 2.3. The operator A generates a Co-semigroup of contractions T(t) 
which is exponentially stable, i.e., there exist two constants K,w > 0 such that 


| Then the 


T(t) |x < Ke", Vt > 0. 


Proof. The first claim follows from [I7] Theorem 7.2.4] and Lemma B.I]and the 
second claim follows from and Lemma [2.2] O 

Remark 2.1. By basic property of Co-semigroups, the operator A*, the dual oper- 
ator of A, also generates a Co-semigroup of contractions T*(t) which is also exponen- 
tially stable. Moreover, the dual system of (Z2) is Y(t) = A*Y(-,t), Y(-,0) € X, or 


yi (a, t) = —Boyo(a, t) + Baya(z, t), 
Yyo(a, t) = — Bry, (x,t), 
ys(a, t) = —Bay4(z, t), 
(2 9) ya(a, t) a — Bsy3(a, t) T Bıyı( a) 
yi(1,t) = aıy2(1,t), y2(0,t) =0, 
y3(1, t) = azya(1, t), ya(0, t) = 0, 
(yi(2, 0), y2(x, 0), ys(x, 0), ya x,0)) = Y (z, 0) 
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2.2. Uniform exponential stability for semi-discrete systems of (2.2). 
In this section we discuss uniform exponential stability for a semi-discrete scheme 
of Timoshenko beam (L-J)-(1.2). The process is the same as [20] [22]. For integer 
N €N, and a partition of (0,1): 0 = zo < a1 < +++ < ay < tN4+1 = 1,2; = ih,i = 
0,1,---,N+1,h = yt. Let vj := (j + 4)h and fi be the value of f(x) at node z. 


— fit fi 


denote the central divided difference operator of first-order derivative f,(v;) and the 
average operator at f(v,;), respectively. Set y:,;(t) = y(a;,t) be grid functions at 
grids {zj}; for 1 = 1,2,3,4 and j = 0,1,---,N. The first four equations of system 
(2.2) are valid at v;, i.e., 


(2.10) { Yr (vj,t) = Boys (vj t) — Baya(vj,t), Yo(vj,t) = Biyi (vit), 
` Y3 (Vj, t) = Baya (vy, t), ya(v;, t) — B3y3(v;, t) + Biyt (vj, t). 


Now replacing y(v;,t) and y;(v;, t) by Uy jtt (t) and bx Yi jth (t) in (2.10) respectively, 
we derive a semi-discretized finite difference scheme of system (2.2) as follows: 


Yaja (Ë) = baðsya j4 (t) — Baya j+ ©) 
Yo.j+4(t) = Pidzyy j+4 (t), 
(2.11) d3,j+4 (t) = baðsYa,j+4 (t), 
Üa jg (t) = 630293544 + Pith +3) j =0,1,: N, 
yı, N+1(t) = —Q1y2,n41(t), Y2,0(t) = 0, 
y3, N+1(t) = —a2y4,n+1(t), ya,o(t) = 0. 


For more information, one refers to 22] and [24]. In the state space C4(N+)), define 


Zir 24,1 Zk,0 
21,2 Zk,1 
(2.12) Xp, = Zh, = Zah € CAN+1) r Zih = i ; Zkh = b , 
Zah : : 
Zan Zi, N+1 Zk, N 


with k = 1,3, l = 2, 4, and the inner product 


4 N 
a (Zn, Znyn =h > (Brzii44% 43) Wa Žn € Xn, 
: l=1 j=0 


22,0 = 24,0 = 0, 21,N41 = —0122,N41, 23,N41 = —Q224,N41- 


To formulate (2.11) into vectorial form, we introduce some matrices. Let An, Bn, 
Ch, and Dp be square matrices of order N + 1: 


1 1 1 -1 


A,=-— i a B,=7 a l 


Ch = diag(0,--- ,0,—a1) and Dp = diag(0,--- ,0,—a@2). Set 
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An 2-1C), 0 0 0 B2B} 0 —p4A} 
gal AF 0 0 vu, | ABa 0 0 0 
a o O Ap a pen 0 0 0 B4B; 
0 0 0 Bi Ap, 0 -BB 0 
0 0 0 
a| I Bih— Ch, 0 
h> Ih o 0 0 
0 2-1BiCh 0 ae 


The state variable of (2.11) is Y,(t = (yın (t) Yan(t ),U3h (t), YAh (t))" with 


yin(t) = (yi,o(t),--- ee yon(t) = (y2,1(t), ++ , ye,n+1(t)), 
ysn(t) = (ys,o(t),--- y3, n (t)), yar(t) = (ya, (t), , ya,n4i(t)). 


Then (2.11) is abstractly written as 


Yn (t) = An Yn (t), 
2.14 
aa { Yp (0) = (yt, y2, y, y)" € Xh, 


where A;, = Dr (Yh + Qa) because ®,, is evidently invertible. 
The classical semi-discrete scheme is similar with (2.14) where the average oper- 
ator P, = diag{ Iyn } 1; IN } 1; IN | 1, IN | 1}, i.e., 


(2.15) { (t) = BrYn (t), 
Yn (0) = (Yh Yk Yn Yh) € Xh, 
with Ba = Yr + Op. 

Here we explain the significance of the discrete scheme (2.14). We plot Figures 1 
and 2, respectively. Figure 1 depicts the maximal real parts of the eigenvalues of the 
classical semi-discrete scheme (2.15) and the semi-discrete p- E e He 
size h. Figure 2 depicts the distributions of the eigenvalues of (2.14) and in 
which N = 100. From Figures 1 and 2 we see that the real parts 7 ses of 
An approach to a negative number and those of By, approach to zero. In both figures, 
we take 8, = b2 = p3 = B4 =a, = aq = 1. 


x without the average operator 
400 x _ with the average operator 


200 
-0.3 100 ed 
ä ppa a ye 
-0.5 “100; 

-200 


“0 50 100 150 200 250 300 350 400 -3 -2.5 2 -1.5 A -0.5 0 
N Real part 


Maximal real parts of eigenvalues 
: 3S : r 
R 
Imaginary part 
o 


(a) Fig.1. Maximal real parts (b) Fig.2. Eigenvalues distribution 


Furthermore, the inner product in X, can be rewritten as 


(Zh, Zryr = h(®,Zn)*Hn©nZn,V Zn; zy, EX, 
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with Hn = diag(81In+1, 62lnw41, 83In +41, baln). The energy Ep(t) of (2.11) is 


4 N 


(2.16) Enlt) = S(Valt), YalO)n = $Y Bley OW 


l=1 j=0 


The following Lemma Pis the discrete counterpart of Lemma [2.2] 
LEMMA 2.4. The operator An defined in (2.14) is dissipative on Xp for all 
h € (0,1). In other words, the solution Y),(t) to (2.11) satisfies 


(2.17) Ey, (t) = —a1 b1 B2|y2,n41|? — 0283Ba|ya,n 41”. 


Proof. It suffices to consider only real solutions. Multiplying the first four equa- 
tions of (2.11) by ABiy,.j42 (l= 1,2,3,4) and adding up for l and j, we obtain 


N N 
Ein (t) = hB1Bo Ds Yaj+400¥2,j43 + do 25+ 3 50 j+ 
j=0 j=0 


N 


N 
+hbapa | > 143,542 5e¥a+4 +> YajthinYs jti 
j=0 j=0 


A simple calculation shows that 


N N 
h |X Yi jpaSe¥ajea +Y Y2 gpa Seti gta | = Y1,N+1Y2,N+1 — Y1,092,0; 
j=0 j=0 


N N 
h 5 Y3,j4+250Y4 544 a 5 Ya,j+4 Ox ¥3,j44 = Y3,N+1Y4,N+1 — Y3,0Y4,0- 
j=0 j=0 
The above two identities, together with y1, N+1 = —a1y2,n+1 and y3, N+1 = —Q1Y4,N+1 


give equality (2.17). O 

Next we show a mixed observability inequality involving the final state observabil- 
ity (see, e.g., Definition 6.1.1]) and the exact observability, which is the discrete 
counterpart of [I7] Lemma 9.1.2] or [33] Lemma II.1] and plays important role in 
verifying the uniform exponential stability of (2.11). 

LEMMA 2.5. Set m = 27! max{1, 6485", Bibs} hq = min{1, mt}, n = 
imh, g = TËT Let M' and M be positive constants and assume that O < h < hı 
and MI4 < H < M'I,. Let the output O(t) for be defined by 


(2.18) O(t) =HYw41(t), Ynsi(t) = (y1, N1 (t), Y2,N+1 (t), Y3, N+ (t), ya, nta (t)) 


Then for any solution Yp (t) to (2.11), there exist a positive constant T such that 


(2.19) flow esdt + CHEKO) > e(h)En(r, 
0 
where the positive functions c(h) and C(h) are given by 
N+1 
n 2 2 
2.2 h) =4mMr— L -L ah H. 
(eae Sn "Cd — 1N+2)(1 4+ mh) 1+mh’ ay) 1+mh 
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Proof. Firstly, for j = 0,1,--- , N+1 we set Y; (t) = (y1,;(¢), y2, (t), ys,j(t), ya (0) ' 


fit) = hy;! (t)HY;(¢) and define the auxiliary functions F}(T) by 


(2.21) B= f hoa 


which is written simply as F}. A simple calculation gives 


(2.22) 
fmi =; f ime- (tat 
= | S allna OP- long (t) Pat = 2h A Yeni Ouija, (Edt. 
l=1 


It follows from (2.11) that 


BiðzY1 j+ (tY j+ Ct) = daj (ur +4 E), 

(2.23) B02 Yo, 1(t)Y2,j i) = Vjt} (t)Y2 j+ a(t) + Baya jg (t)Yo542 E), 
P3ðzY3 j+4 (t)Y3 j+ (t) = a j44 uss a(t) — B1Y1, j3 (Cys 544 (6), 
BabxYa,j OE a(t) = 93.544 (t)Ya j+ ©), 


which means that 


T4 
| 5 BiðsYi j+ (Yi j+ (tdt 
l=1 


d 
= | Fo Owy O + vaga Om seg Ole 


2.24 T 
a z f Bayn j44 (Yara (Ë) — Bada gya (EUa gya (Cat 


T 


= ERETO) + Ys j+ Ya jti (t) 7 


= 
-f Pry j+ (tys, j+ (t) = Baya ja (tua, j+ (t)dt. 
0 


Combining (2.22) and (2.24), we arrive at 


FaF; r 
+ =2h Eee (t)Yo 542 (t) + Y3,j+4 (t)Ya,j+2 (t) g 
(2.25) 


-on f Bit ja Os j+4 E) = Baya gra (O)Yo,j44 (t)dt. 


Secondly, for the second term of the right-hand side of (2.25), we apply the 
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Cauchy-Schwartz inequality to obtain 
—2h L B1Y1 544 (t)Ys,544(t) — Paya j (tY, j+ (tdt 
> =l [ Br (lng Ol? + lus Ol? + lyr z+ Ol? + lys541)/7] dt 
-} L Ba [lyzi ©? + lya OI? + luz jt)? + lya j+ (t)|?] dt 


h T 
=F f eulmOP + FEBalans (0)? + 52 Bal + Balun s(OPae 


h [7 B E 
| balyi jat)? + Æ bolua, j1 (6)? + Eblus z416)? + Balyages (t) Pat 
0 


Bo Bs 
2 —m(Fy41 + F;). 


Substitution of it into (2.25), for j = 0,1,--- , N, we obtain 


(226) Fan > Fj + 2h [vn ja Evoje + vaga Ouar O 
where 7 = i Tmk h and 6 = 7 TËT By deduction argument, we obtain 
Fj(r) < WD Fyi (T) 
N T 
(2.27) = 2nd X nF) Eee (tyz j+ (t) + Ys jra (tya ja E) = 


k=j 


Finally, since from Lemma[2.4] En (te) < En (tı) for any tı < te, we deduce that 


L En (t)dt > En(T) f dt = TE, (rT) 
0 0 


This, invoking definition of F; and E;,(t), and the estimate (2.27) leads to 


r N+1 


2rEx(r) <2 | moa f S S viet (8) Pas f Y fy (tat 
0 9 j=01=1 j=0 
N+1 —(N+1) _ 
-Y ys tain 
; =n 
j=0 
NHIN s 
-2h XO on ln jeg Our 13) + vag Oua O, 
j=0 k=j 
g+) — N+1 es : 
= 1-7 4 Py salt ao Dr 7 Dust yP + |yi,n4.2 (0)| 
nN +D) = 


< 1-7 1 [Fn (T) + 26M~'(En(r) + En(0))] , 
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which implies that 2r—wsh— En(r) — 29M~\(En(7) + En(0)) < F4i(r). This, 
together with the following inequality Fy4i(r) < hM! fy ||HYn41(8)||?dt, gives 


(2.19). O 
Remark 2.2. By the definition of n and (N +1)h = 1 in Lemma[2.3] it is easy to 


see that 
li _ ii 1—mh 1—mh T ói 
im = lim =e ; 
h—0 n” h>0\1+ mh 1+mh 
This means that 4mMt aon is uniformly larger than aR for all h € 


(0,h2) with some hg < 1 whenever we choose T large enough. Moreover, since 
Tok < 4, we can choose T to ensure that c(h) and C(h) defined by (2.20) satis- 


fying Supp €(0,h>) on <1. 

Now we are in a position to state the main result of this section. 

THEOREM 2.6. Let h* = min{hı, h2} where ho is specified in Remark B.A and 
hı by Lemma [B5], and let T;,(t) be the semigroup generated by An defined by (2.14). 
Then T;,(t) is uniformly exponentially stable with respect to h € (0,h*), i.e., there are 
two constants K and w independent of h such that 


(2.28) IITn(4)¥n(O)|Ix, < Ke Yn (Oln V¥n(0) € Xn. 


Proof. By RID, yi,w+i(t) = —a1y2,N+1 (t) and y3,n4i(t) = —a2y4,N+1(¢), it 
follows that 
(2.29) 


Er(T) D= f Én En(t a=- f 0181 B2|yo,n+1(t)|? + a283b4lya,n+1(t)| dt 


Sei pe HY 41 (Ċ)||?dt, 
0 


where k = min [tt Sot. sth: a } . Plugging (2.19) into (2.29), we obtain 


En(7T) — En(0) < —k[c(h) Ep (7) — C(h) Ep (0)]. 
Thus E(t) < THON En (0) and so SUPhe(0,h*) AA < 1 for some T by virtue 
of Remark 2.2} By Co-semigroup property, this turns out that the energy Ep(t) of 
system (2.11) decays uniformly exponentially to zero as t + oo. In terms of the 
relationship between ||Y;(t)||x, = ||Th(£)Yn(0)||x, and energy E),(t), we have finally 


obtained (2.28). O 


Remark 2.3. Similar to Remark[2.1| the dual system 


Yn (t) = Až Yn (t) 
2.30 nant); 
oe { y(0) = (Yn, y2 UR YA)! E Xn, 


which is 
tat) = —B2d2Yo gay Paya Jj+4 (t), 
Yoj+4(t) = —Brdeyy grat), 
Y3,j+2 (t) = —646 zya grt), 
(2.31) Ua j+ (t) = —Badoya j+ (t) + Fir j+ (t), 9 =0,1,---,N, 
yı, N+1(t) = a1ye,n4i(t), Y2,0(t) = 0, 
y3,n+i(t) = azya,N+1 (t), Y4,0(t) = 0, 
Yn (0) = (Yh: Yn: Yh Yn) E Xh, 
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is also uniformly exponentially stable. 


3. Convergence of solutions of to (2.2). The objective of this sub- 
section is to show that the Co-semigroups 7),(t) generated by Ap defined by 
approximate T(t) generated by A defined by in some sense. 

THEOREM 3.1. (Trotter-Kato [16)) For every h € (0,1), assume that there 
exist bounded linear operators Pn : X — Xn and Fh : Xn > X satisfying 

(Ai) There exist two positive constants Mı and Mə such that ||F),|| < Mı and 
| Pall < Ma; 

(A2) | Fr Pay = Y||x >O0ash—>0 for allY € X; 

(A3) BE = In. 

If (A1) and (A3) are fulfilled, then the following two statements are equivalent: 

(a) There exists a Ap € p(A) NMP, p(An) such that for all Y € X, 


|| Fh Qoln — An) PAY — (Aol — A)~'Y||x 7 0 as h > 0. 
(b) For every Y € X and t > 0, 
(3.1) ||Fn Tp (t) PhY = T(t)Y ||x >O0ash>0 


uniformly on bounded time t-intervals. 

In order to apply Theorem [B.I] one always replaces property (a) by a condition 
involving convergence of the operators Ap to A in some sense [16]. 

PROPOSITION 3.2. Let the assumptions of Theorem [3.1] be fulfilled. Then s- 
tatement (a) of Theorem [3.1] is equivalent to assumption (A2) of Theorem [3.1] with 
additionally the following two conditions: 

(C1) There exists a subset D C D(A) such that D = X and (AoI — A)“ D = X 
for ao > w; 

(C2) For all Y € D, there exists a sequence Yp € D(An) such that 


(3.2) lim FrYp, = Y and lim F, Any pn = AY. 
h-0 h-0 


To apply Theorem [B.I] and Proposition B.2] to the convergence analysis in this 
section, we firstly introduce the operators Fn : Xa + X and Pr : X 4 Xp. Let xs be 
the characteristic function of a set S. Then Fp is defined by 


N 
pane Yrit4x 


i (xi,vita] Yih 
FY, = ico Y2it+bX(eivi41] WY, = Y2h EX, 
(3.3) Da Y3 i+ 4X (21,2; 1] Y3h 
N 
Vigo Ya itd X (ei ti] Hah 
Y2,0 = yao = 0, Y1, N+1 = —Q1Y2,N+1; Y3, N+1 = —Q2Y4,N+1; 


in which the second row was assigned to unify the notation of y; j4 L On the other 
hand, P, is defined as, Y Z = (z1 (£), z2(a), 23(x), za(x))' € X, 


(3.4) P,Z = (T'a (2), I?z9(2), z3(2), I*za(a)) ' , 
where 
Hals) IF zx(2) 
zale) Ik zy (x) 
Pala = 2h] 7. , T¥z,(x) =2h-1| 7 (418, k=2,4, 


I zala) Tr qe) 


202301.00193v1 


chinaXiv 


CONTROL CONVERGENCE FOR A TIMOSHENKO BEAM 13 


and 
EN+1 
I} z (x =f a)dx + ay I} 4122(2), 
Tj+1 
Iz (z F — T}4121(2) (j = 0, 1, a, N- 1), 
I? zo(x) = [al z)dx, I?zo(x) = / zo(x)dax — I?_,2z0(x) (i = 2,3,---,N +1), 
Lia 
@N4+1 
I, z3(x =f x)dx + ag ly, 1 24(x x), 
DP. 
I3z3(x ya fe x)dx — I 3123x) (j =O, 1, <, N-1), 


Page- m aa, Pa (eT ey BHI NSD, 


LEMMA 3.3. The operators F, defined by and P, defined by satisfy 
assumptions (A,)-(A3) in Theorem[3_]] with Mı = Mz =1. 
Proof. By the definitions of the inner products of X and Xp, it is easy to see that 


2 2 
(3.5) | FnYnl? = Dy! cB Insa) = Yl, 
which implies that ||F}|| = 1. Moreover, by setting Jf = Ij = 0, Ij, = —aily44 


and ETS = -alk 44; we have, 


4 N ə N iin 
IZR, =a HY e| af laea < [ZR 


l=1 i=0 l=1 i=0 Ti 


This shows that ||P,|| < 1 and hence assumption (A1) holds. 

To prove assumption (A2), set D = [C1 (0, 1)]*. It suffices to show that || Fn PRY — 
Y||x > 0 for all Y € D since F, and P, are bounded linear operators and D is dense 
in X. In fact, for Z € D we have 


i+1 i+1 
| Fa [Pr Z] — TEDD [ lzi(r) — a(x)? drdz. 


By the same idea of [16] Section 4, p.34], we obtain that (Az) holds also. Finally, 
using Y2,0 = y4,o = 0, we have 


N 
Yk,j é 
(3.6) I} pres cee = h k = 2,4, J= 1,2, ete N +1. 
i=0 


Substituting (8.6) into Z! a Wit t Xens] (l = 1,3), we obtain 


N 
Yl, j P 
(3.7) ij ye iene = n=, l= 1,3, J= 0, 1, 7 N, 
i=0 
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where y1, N+1 = —Q1Y2,n+1 and y3, N+1 = —Q2y4,n+41 were used. (8.6) and (8.7) then 
imply (A3). O 

THEOREM 3.4. Let T(t) and Th(t) be the Co-semigroups generated by the opera- 
tors A and An which are defined by and (2.14), respectively. Let Y(0) € X be 
the initial value of 2.2) and set YP := P,Y(0) be the initial data of (ZII). Then 
F Yp converge to the initial value Y(0) in the sense of Fh Yp > Y(0) as h > 0, and 
|| Fn Tp (t)¥;2 — T(t)Y (0)||x 4 0 as h + 0 uniformly on bounded time t—intervals. 


Proof. The condition (C1) of PropositionB.2]is obviously valid whenever we choose 
D = DA D(A) and ào € C with Redo > 0. It is therefore sufficient to prove (C2) 
of Proposition B.2] Indeed, for Y = (yı (£), yo(x), y3(x), ya(x))' € D, we introduce 
Yn € Xp through 


Tin = (yi(@o), yi(ti), «++, yi(@n))', Fon = (Y2(21), y2(w1), «++, yo(en41))", 


Tan = (ys(z0), ys(a1), =, ys(tn))', Gan = (ya(w1), yi(z2), +++, ya(en))", 


and Yr = (Yin, Y2n: Y3h Van) | We thus have 


N Ti+ N x 2 
Zizo i yilt ATX (wi wipı = Žo 2,4 T)ATX (x; 2:41] 
N Tit ae 
||FrYn _ y|? _l Vigo Ja ¥2( T)ATX(x; 15 
= 4 2a SP u(r OPN pees] 
4( MPN iy esa b 


N it N 
Jio J: “yf (T ATX (a; i41] — Da0 je Ya\T 


N z 
5 Í Y(T)ETX (4,041) 
i=0 “Ti 


2 


N 2 


2 
+ dx 


IA 


Ti+. 
1 
/ Yı Cae cements 


Ti+. 
J yi(r)dr 


ie 1 
Def i=0 
1 4 2i 
=>) 8 

IDAJ È 


2 g 
X(£i,£i+1] F / yı(r)dr 


Kere) dx 


2 


1 4 N Dit. Lit. 2 1 4 N Lit. x : 
1 
= DaD dx i. y(t )dr +a dx I yilT)dr 
i= i=0 Y Ti l=1 =0 “Ti i 
By the Cauchy-Schwartz inequality and Fubini Theorem, we obtain the following 
estimates 
N Lit. Lit. 2 N Lit. Lit. 
S| af oal E [eis ade [ery ar 
i=0 Ti z i=0 " Ti x 
2 re E 1 2 2 a oa 1 2 PATI] 
=Y af em-oa S MOP dr = PON 
i=0 7 Ti Ti i=0 7 Ti 


i x 2 . 
Similarly, we have par o dx} f yi(r)dr < k’?llyi(T)llz2,} = 1,2,3,4. Combin- 
ing three inequalities above, we obtain Fa Yn > Y as h > 0. 
Set Zn = AnYn € Xn and we obtain ©), Zn = (Wp, + Qa)Y n. By the expressions 
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of ®,, Yh, and Qa, we have 


Litt) — Y2rvi Ti + ya(x; 
#iged = py 22! 2- y2(zi) — p Ean tut, 
r= yi (xi41) — Yı (xi) 
2.j+3 $ 
(3.8) Tra h 
_ a YalTi+1) — Ya(zi) 
23.542 = Ba ; l 
Z4 gia E ys (xi) + Tan talea, j=0,1, N, 
where we set artificially that 22,0 = 24,0 = 0, 21,v41 = —@122,w41 and 23,v41 = 
—a22z4,N+1 to unify the notations Zl j+} for j = 0,1, --- , Nn and l = 1,2,3,4. Plug- 


ging into the operator F, and the same techniques as above, it is easy show that 
Fa Zn —> AY as h - 0. For the consideration of the page limit, we omit the details. O 

Remark 3.1. For the Co-semigroups T} (t) and T*(t) given in Remarks [2.1] and 
we have the convergence 


|F TOY — T*(t)Y (0)||x > 0 as h > 0, 


uniformly on bounded time t-intervals, where Y(0) € X is the initial value of 
and YP = PhY (0) the initial data of : 


4. Convergence of the controls. In this section, we give exact controllability 
of (I), in particular the analytic form of the control that drives the initial state to 
rest. Same to (2.1), we introduce the transforms 


a(x, t) a Wa(x,t) = olz, t), q2(x, t) = pwr(x,t), 
q3(x, t) = Qalx, t), qa(a, t) = pot (2, t), 


and obtain equivalent form of (1.1) 


d(x, t) = B2q3(x, t) — Baga(z, t), 
g(x, t) = b1qi (x,t), 
G3 (a, t) = Baga (a, t), 

(4.1) qalx, t) = B3q5(a,t) + Big (x,t), 
q2(0,t) = 0, q4(0,t) = 0, 
qı(1,t) = u(t), q3(1,t) = u2(t), 
a(x, 0) = q! (x), l = 1,2,3,4, 


where 6; (l = 1,2,3,4) are explained after (2.2) and q! (£) = w(x), ¢?(x) = w(x), 
(x) = $°(x) and q*(x) = ¢'(x). The aim of this section is twofold. The first one is to 
design controls u(t) and u2(t) which drive the solution of (£1) from the initial state 
to rest at finite time 7. The second one is to construct a family of discrete controls 
which converge to the continuous one. We split this section into two subsections. 


4.1. Exact controllability of (4.1). The state space of (4.1) is X also. Set a1 
and az to be zero in Wg defined by (2.6): 


0 BB 0 0 BEO 0 0 
1 —1 0 0 0 0 1 0 0 

4.2 Wg = — z is 
mi. ammo 0 0 fbx! 0 vA 0 
0 0 —1 0 0 0 0 1 
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The corresponding operator A defined by (2:5) is denoted by Ao in this case: 
AoY (x) = Pi (HY (x))' + Po(HY (x)), Y Y(x) € D(Ao), 
4.3 
ne) D(Aq) = ÍY € X: HY € H4(0,1}:C%), we | =o}, 
It follows from Lemma[2.2|that the operator Ag is skew adjoint. Now we consider 
the dual system of (41), which is corresponding to (— Apo, B*), i.e., 


P(z,t) = —AgP(z, t), P(z,t) = (pı (x,t), polz, t), pa (x, t), pa(a,t))', 
wa) d PEO = P'a) = (P(e) pe) 282), vi)” €X, 


* 0, 81626(s — 1), ,0 
O;(t) = B* P(x, t) = Gare S a P(a,t), 
or 
Pilz, t) = —Bop9(a,t) + Bapa(z,t), 
po(x,t) = —B1p) (x,t), 
p3(x,t) = —B4p4(a, t), 
(4.5) pa(x,t) = —B3p3(x,t) — Bıpı (x,t), 
pi(1,t) = p2(0,t) = ps(1, t) = pa(0, t) = 0, 
pi(x,0) = p'(x),1 = 1,2, 3,4, 


O71 (t) = (61 Bepe(1, t), B3Bapa(1, ae 


where O(t) is the output of adjoint system of (4.1) without controls, 6(- — 1) is the 
Dirac delta function, and B* € L(X,, C?) is the observation operator. The boundary 
control system (4.1) can be reformulated as follows 


s Q(-,t) =A Q(t) + Bu(t), 
eee Naas { Qle,0) = O(a) = (a'a), (2), q(x), ¢4(a))7 € X, 


where Q(x,t) = (q1(a,t), go(2, t), q3(x, t), qa(a,t))' is the state of system, u(t) = 
(ui(t), u2(t))' € L2,.([0,00),C?) is the input and B € L(C?,X_j) is the control 
operator. To apply Theorem 2 of [II] to obtain the exact controllability of (ZI) or 
(£6), we first present an admissibility result. 

THEOREM 4.1. Let To(t) be the unitary Co-group generated by Ao defined by 
(4.3). Then the control operator B defined in is admissible for To(t). 

Proof. It suffices to consider the real part of the solution. In light of duality, 
we only need to show that the observation operator B* is admissible for the semi- 
group Tj (t) since (Ap, B) and (— Ao, B*) are a pair of dual systems. Now, for every 
P°(x) € D(Ao) and the state trajectory P(x, t) = Tj (t)P° (x) of (£4), we introduce 
the auxiliary function 


Pit) -j xpı(x,t)polx, t) + xp3(a,t)pa(a, t)dz. 


Differentiating P(t) gives 
(4.7) 
P(t) 


= f [xpi (ax, t)po(x, t) + wpi(a,t)po(a, t) + xps (x, t)pa(x, t) + xp3(a,t)pa(a, t)|dx 


4 


1 
= -f 2 trp; (x, t)pi(x, t) — Baxpo(x, t)pa(x, t) + Bıxpı (x, t)pa(x,t)| da. 
0 


l=1 
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Performing the integration by parts gives 


4 1 
49 f Bwile plede = Zll + alpa AP) — BoC), 


and 
1 
(4.9) | Baxp2(x, t)pa(a, t) ~~ Bi xpi (x, t)ps(a, t)dx < mıEp(t), 
0 
in which 
il 
(4.10) Ep(t) = 5llPC ÐI 


is the energy of (4.4) and mı = max | yB Aa, Vara}. Combining equations 
(£7)- (£9), we derive 


(411) PA < -ibl t)? + Balpal1, DP] + (1 +m) Br). 
Similarly, 
(4.12) |P(t)| < m2Ep(t), Vt > 0, 


with mo = max | Ar r, Verte" . Integrating (£11) from 0 to 7 and using 


(4.12) yields 
(4.13) 


[Pma + balpa(1, t)|?]dt < af (1+ m1)Ep(t)dt + 2m2|Ep(T) + Ep(0)]. 


Finally, using the expression of Oı(t) and the conservation of the energy Ep(t), 
we obtain 


(4.14) T |O1(t)|?dt < k-Ep(0) 


with k, = 22Gtm)+?m2) and m3 = min {8,'8,*, 87'863’ }-This means that B* is 
admissible for Të (t). O 

In light of the infinite dimensional version of Russell’s “controllability via stabil- 
ity” principle, we get the control which steers the initial state of (46) or (£1) to rest 
in some time T. 

THEOREM 4.2. Let T(t) be the Co-semigroup generated by Aj — BB* and the 
time T is chosen so that ||T(T)||x <1. A control u(-) € C([0, 7], C°) for (£9) or F) 
driving Q°(-) € X to rest in time T is given by 


(4.15) u(t) = —B*[W;(-,t) + Wo, t)], 


where Wy = (w1, we, w3, wa)! and Wp = (wor, Wo,2, Wo,3, Wo,4)| are the solutions of 
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a forward system 


wi (2, t) = B2w3(x, t) = Bawa(2, t), 
w(x, t) = Biwi} (x,t), 
w3(x, t) = Bawy (x,t), 
(4.16) wa(x,t) = 83w4(x,t) + Biwi (x, t), 
w2(0,t) = 0, wa(0,t) = 0, 
wz1(1, t) = — BG 2we(1,t), ws(1, t) = —ß3ß4w4(1,t), 
w(x, 0) = w (£), l = 1,2,3,4, 


thy,1 (x, t) = Bow, a(x, t) — Bawe,a(z, t), 
wp,2(2, t) E Piw, (a, t), 
Wb,3 (x, t) = Bawy a(z, t), 
(4.17) wo,a(x,t) = Bzw, 3(2, t) + 61w 1 (x,t), 
wb 2(0, t) = 0, wpb 4(0, t) =0, 
wp 1(1,t) = b1b2w,2(1,t), w, 3(1,t) = 8384W, 4(1, t), 
2 


respectively, with W;(x,0) = WP (x) = (I — L-)*Q°(x) and L, = T*(T)T (T). 

Proof. The proof was actually given in proposition 2.2 of and theorem 2 of 
[II] which relies heavily on the admissibility result Theorem [4.]] and the exponential 
stability of the semigroup 7 (t). Here we just give a sketch of the proof for complete- 
ness. 

Firstly, we notice that the operator Ag — BB* corresponds to the operator A 
with a, = 8182 and ag = $384. Furthermore, Theorem [2.3] tells us that T(t) is 
exponentially stable. Thus system (4.16) can be written as an abstract evolution 
equation 


{ Wiet) = (Ao a BB*)W;(-,t), 
W;(z,0) = Wẹ (2). 


Similarly, (4.17) can be written as an abstract evolutionary equation 


{ W(t) = (Ao + BB*)W,(-,2), 
W(x, T) = Wy(a,7). 


Let Q(x, t) = Wy (x,t) — Wi (x,t) for t € [0,7]. Clearly, Q(x, t) and u(t) satisfies (ZI). 

Secondly, let Wp (x,t) = T(t)W?(x) and Wi(z,t) = S(t)W)(a) with S(t) being 
the Co-semigroup corresponding to (4.17) and WẸ (x) the initial value. It then follows 
from (416)-(Z1I7) and W, (z, T) = Wp(x,T) that Wi(a,t) = S(r — t)W, (2,7) and 
OO) = W? (a) — Wẹ (x) = W? (a) — S(T)W; (x, T) = W? (a) — S(T)T(T)W? (2). The 
proof will be accomplished if we can show that S(t) = T*(t). Indeed, on the one 
hand, it has 


(4.18) Wilz, t) = Ax S()W? (2) 


with the infinitesimal generator Ax of the semigroup S(t). On the other hand, from 


(£17), we have 


W,(a,t) = —(Ao + BB*)S(7 — t)W, (2,7) 


(4.19) = —(Ao + BB*)S(t)W2(x) = (Ao — BB*)*S(t)WP (2). 
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Identities (4.18) and (4.19) imply that Ax = (Apo — BB*)* and hence S(t) = T*(t). 

Since 7 (t) is exponentially stable and so is for 7*(t) (Remark B.I). Since there exists 

T such that ||7(7)|| = ||7*(7)|| < 1 and therefore the operator I — L, is invertible. 
Finally, since W, (x, 7T) = Wp (z, T), Q(z,7) = 0 is trivially valid. O 


4.2. Uniform controllability and convergence of discrete controls. We 
begin with (41) and discretize it using the same discrete scheme as that of subsection 
We therefore obtain a family of finite-dimensional control systems: 


Gi j+4 (t) T P2ôxq2 j4+4 (t) — Bada j+ (Ë), 
do j43 (0) = Biðst jt) J=0,1,: N, 
ås, j+4 (t) = Baĝsda j3 (t), 

(4.20) qa, galt) a $ 36093542 (é) + Pid 544 (t), 


q2,0(t) = 0, qa,o(t) = 
M,N+1(t) = uin (t), 93,n-41(t) = uzn (t), 
q = q (ih), L= 1,2,3,4, 1=0,1,---,N+1. 


l 
D 


The energy of system (4.20) is defined by 


N 4 
TODD (lar j+ (t )12). 


g=0 {= 


(4.21) Ec,n(t 


wie 


It follows from Lemma[2.4] that Ec.n(t) satisfies 
Eon(t) = Bi Bourn (t)q2,n-+41(t) + 83 Bauan (t)¢a,n 41(t). 
This motivates us to introduce the output for as follows: 
(4.22) On(t) = BrQn(t) = (618292, +1(t), 8364ga, N41 (¢)) |, 
where the unknown quantities are Q(t) = (qın (t), q2n (t), gan(t), Gan(t)) | with 


qın(t) = (q1,0(t),--+ ,41,w(t)),  Gan(t) = (2,1 (t), , 2, 41(t) 
asan(t) = (g3,0(t),--: ,¢3,w(t)), dan (t) = (qa (t), +, ga,w4i(4)). 


This makes system (4.20) be passive, i.e., Ec a(t) = Ofun (t) and u(t) = (win (t), uzn (t))". 
Theorem [2.6] tells us that the resulting closed-loop system under the negative propor- 
tional feedback u, (t) = —Op(t) is uniformly exponentially stable. Using these basic 
facts, we construct controls up,(t) which drive the solutions of (4.1) from the initial 
data to rest in some time 7. We still apply Russell’s “controllability via stability” 
principle. Now, consider the forward system 


Wi jth (t) = B2bxWe 542 i(t) — Baw, 544 (t), 
tw 542 (t) = Pi d,W, jet), 
ts 54a (t) = Bað Wwa 542 (t), 
(4.23) tg 5428) = abn 542 (t) + 1w j4 (t), j =0,1, N, 
wi,n+i(t) = -ow N+1(t), w2 o(t) = 0, 
w3, N+1(t) = —a2qwa,n4i(t), wa,o(t) = 0, 
win(0) = (wir, ,wiw4i)' E€ CNH, 1 = 2,4, 


wen (0) = (wio, wpn)! E CNH, k =1,3, 
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and the backward system 


tyr j (t) = Boða Wya j+ (t) — Bawa 542 (Ë), 
tya j+ (t) = Bið Wir j+ (t), 
W3 j+} (t) = Badz Wo4 j+ (t), 
(4.24) tpa jz (t) = B3ðz W3 j+4 (t) T Biwi, j+ (t), j =0,1, LN, 
wo1,n+1(t) = a1Wo2,n41(t), wo2,0(t) = 0, 
wy3,N+1(t) = a2wia, N+1 (t), weao(t) = 0, 
wen (T) = win(T), l = l; 2 3,4, 


where aı = 612 and az = b384. The unknown Wp and Wy», of (4.23) and (4.24) 
are defined similarly as Qn, respectively. Hence the desired control u;,(t) is given by 


(4.25) un(t) = —By[Wen(t) + Won (t), 


provided that we choose the initial value of Wta (0) as Wen (0) = (In — L7.n)7'Qn(0) 
with ||7;*(7)||x, < 1 and ||Z;,(7)|lx, < 1 holding for h € [0,h*]. Here T),(t) is the 
Co-semigroup defined in Theorem 2.6] For simplicity, we use the notations Lr n := 
T;(r)Ta(r) and L, = T*(r)T (T). 

Now we present main result of this section, i.e., the discrete controls up,(t) given 
by (4.25) are convergent to u(t) constructed by if some convergence conditions 
are satisfied for the initial data. 

THEOREM 4.3. Choose T so that ||Lr alix, < 1/2 and ||L,||x < 1/2 are satisfied 
for all h € [0,h*]. Let Q? € X be an initial value of (Z-6) and set Q9, := P Q? to be 
the initial data of (4.20). The following results hold true: 

(i). We, given by We, = L- nQÌ (0) converge to WP = £L-Q® in the sense of 


|F Wg, — W2(0)\|x > 0 as h > 0. 


(i). Wer given by Wer = (In — Lra) tQha(0) are convergent to W® = (I — 
L,)~*Q® in the sense of 


| Fn W?n —W°(0)|lx > 0 as h > 0. 


(iii). The discrete controls up(t) given by (4.25) corresponding to W?,, are con- 
vergent to u(t) constructed by (4.15) related to W® in L?([0, 7], C4). 

Proof. By Theorem[3.4Jand RemarkB.1] for arbitrary Q € X, as long as || F, PLQ- 
Q||x —> 0(h > 0), there hold 


(4.26) ||FnTh(7)PrQ — T(7)Q\|x > 0 and || Fi, TF (7)PrQ — T*(7)Q||x > 0. 


Let u = SUPhejo,n] (TE) Zn (7)Ilx,)- Then y < 1 because ||T;(r)||x,, < 1 and 
\ITn(7)Ilx, < 1 hold for h € [0, h*]. 
(i). Since Pa Fh = In, by the continuity of T*(t), 
(4.27) 
Fp WO, — W2(0) = Fab nQ) — LQ? = FuT#(r)PalFaTa(7)PaQ?] — £Q? 


= FT} (r)Pr[FrTn(7) P,Q?) — T* (7) Fa Th (7) Pn Q?| + T*(7)[FaTh (7) PnQ®*| — £,-Q°. 


By (4.26), both the sums of the first two terms and the last two terms of (£27) are 
convergent to zero as h — 0. This proves (i). 
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(ii). Firstly, for positive integer n, it is easy to obtain that 
FL? nOn — LEQ° = Fab? Pal En Len PhO] — £21 [Filer PhO] 
(4.28) + LI HFa Lr n Ph] — LITH’). 
By induction, u < 1 and the identity (4.28), we know that 
(4.29) FRL? pn — L2Q°|Ix > 0 as h > 0 Yn EN. 


Secondly, we prove that the convergence of (4.29) is uniformly for n € N. Indeed, 
by assertion (i) just proved, for € > 0, there exists h’ € (0,h*) such that 


| Fn Lr nP? —£-Q° lx < - for all h € (0, h’). 
Thus, by Theorem [4.3] since ||7*(7)||x < 1 and ||T(7)||x < 1, it has 
E 
(4.30) [L7 EhLrn Pr? — L-Q°]llix < E xl ELP? — £-Q" IIx < 5 


for every positive integer n. On the other hand, since ||Z-,,||x, < 1/2 and ||£-||x < 
1/2, we can select a positive integer k such that 


E 
ToT E< ao" 


which holds for all n > k. It then follows from (4.28)-(4.31) that for n > k and 
h € (0, h’), 


(4.31) FLL; Pallx < ——— 


Th 


|| Fn Le Qh — LEQ IIx < (FAL Pall + IERT x) Fa LenPrQ IIx 
+ ee r,hPaQ? — a QIIE IIx 


(4.32) aeol’ + 5 


TOT 2 


Choose h” € (0,h*) so that 
(4.33) |FL nAn — LQ IIx <e,n=0,1,---,k-1,he (0,h"). 


(432) and (4.33) imply that Fp L? rh ? are uniformly convergent to L? Q° for all n € N 
as h > 0. 
Finally, by basic theory of analysis, we have that 
lim [Fp Wh — W° (0)] = lim 1 [Fn (In — Len) PRQ? — (I — £)'Q*| 
= 
= lim X [Fh L} p PQ? — L2Q => lm lim [Fh (Le h)” PQ? — £2Q°| = 


h-0 


which show that the second assertion (ii) holds. 
(iii). Consider the mixed variable 


(4.34) Mn (t) = FrlWen(t) + Won(t)] — [We (t) + Wi (t)] 


in which W(t) and W,(t) satisfy (4.16) and (4.17), and Wyp,(t) and Wpn (t) satisfy 
(4.23) and (4.24), respectively. Thus Mp(t) satisfies 


(4.35) { Mn (t) = A[FnW pn (t) — W(t)] + A* [Fn Won (t) — Wo()], 
l Mn (0) = Fa [W pn (0) + Won (0)] — [W (0) + Wa (0)]. 
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Let Em, (t) be the energy of Ma (t). Then 


Fult -5f MT OHM, (t)at = MON: 


It follows from the same calculation as that of Lemma B.2] that 


Em, (©) =—B{B3|w2,n41(t) + wez, ny1 (t) — w2(1, t) — w,2(1, t)? 
— 6383 |wa,n+i(t) + wea,w4i(t) — wal, t) — we,a(1, t)’. 


By the expressions of u, (t) and u(t), we obtain 


[iaw )—u (t)||Z2dt 


= pib? [ [w2 w4i(t) + wee,w41(t) — we(1, t) — woo(1, t)|?dt 
tf |wa,nw41(t) + woa, Nyilt) — wa(1, t) — wo,a(1, t)|?dt 


~~ [Bra (t)dt = Em, (0) — Em, (7) 
1 


1 
= sllFaWin — W° + Fa Win — We ll - 5llFnTh(7) Wp, —T(r)w® 
+F Ty (T) Win — T* (7) We Mlk. 


By the assumptions and assertions, we finally arrive at ||un(t)— u(t)||z2 > 0 ash —> 0. 
Oo 

The inverses of operators I» — L+,» and I — £, are involved in Theorem [4.3] and 
are given by Neumann series. It is well known that the curse of dimensionality may 
happen. In this case, the numerical approximations of controls given in Theorem 
[4.3] cannot be realized from computational point of view. To get off this hook, we 
truncate Neumann series of the inverse of In — L,,, and I — £- in light of a prior error 
and further obtain precise error estimates between the discrete system controls and 
continuous system ones. We elaborate this problem in the following Theorem [4.4] 

THEOREM 4.4. Let 7, ra E€ X, Qn = PrQ°, Win: and Wẹ, be the same as 
those in Theorem cA Set W, hes to be Yo" o[Lrn]"Q?. Replacing the initial value of 
(4.23) 4.23) by W, he n: the corresponding solution of (4.23 (4.23) is denoted by Wn h(t). Then for 
any € > 0, there exist hn € (0,1) and n € N such that the approximating controls 
Un,r(t) = —Bs[Wn,n(t) + Wy,n(t)] satisfy 


(4.36) L banawi e for all h € (0, hn). 


Proof. It is easy to see that 


(4.37) pi [tin a(t) -ulat <2 f jatt lèadti+2 f lun (t)-u(t)l|22dt. 


Based on assertion (iii) of Theorem [4.3] there exist hn € (0, h*) such that 


(4.38) f lur (t) — u(t)||2edt < - for all h € (0, hn). 
0 
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Furthermore, by definitions of un,a(t) and up(t), we consider variables 
(4.39) Mn, n(t) = [Wn a(t) + Won(t)] — [Wn (t) + Won (t)] 


which satisfy 


Mn x(t) = An Mn a(t), 
(40) ae .(0)— Wo, W%,. 


On the one hand, we have 


[Mnallx, = |T OW, n -Wallig = IT 


1) SD Ling 


i=n+1 


Xp 
< Ky (1-9) QR 


where K is given in Theorem [2.6]and 7 = || L+ al|x,. This means that when 


-1,, vE- n) 
> (lnn) t In > 1, 
2K Qh llx, 
there holds 
vE r 
(4.41) [Mn a Olx, < F for any t € [0,7] and h € (0, h*), 


On the other hand, if let Em, „(t) be the energy of (4.40), then we have 


1 1 
Em, n (t) = 5 (Mn, (t), Ma (Ox, = 5ll Man Olli, 


From the same calculations as that in the proof of assertion (iii) of Theorem [4.3] we 
obtain 


Em,,,, (t) = —8783|Wn2,n41(t) — we, +1(t)|? — 387 |Wna,wai(t) — wany (t). 


By the expressions of un,,(t) and u,(t) and using (4.41), we finally obtain 


(442) f |juna(t)—un(t)leedt=— f Én (Edt = Eas, (0) ~ Bat, a(7) < É. 
0 0 


Combining (4.37)- (4.38) and (4.42) gives (4.36). O 


5. Concluding remarks. In this paper, we propose a general finite-difference 
based semi-discrete scheme for a Timoshenko beam under boundary control modelled 
as a port-Hamiltonian system. A key feature of this semi-discrete scheme is that 
the discretized systems are still port-Hamiltonian systems and possess some common 
properties of the continuous system uniformly. This advantage is reflected from dif- 
ferent perspectives. Firstly, it allows preservation of certain properties such as energy 
conservation and control design from infinite-dimensional model to finite-dimensional 
approximation. Secondly, this discretization preserves the uniform exponential sta- 
bility with respect to the step size. Thirdly, the proof for the uniform exponential 
stability of the semi-discrete systems is parallel to the continuous counterpart, which 
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is somehow an indirect Lyapunov method and is different to existing approaches be- 
cause the Lyapunov is not directly positive definite and must be combined with the 
Co-semigroup property. More importantly, the semi-discrete systems are uniformly 
exactly controllable and the discrete controls are strongly convergent to the contin- 
uous counterpart. Although we focus only on the Timoshenko beam equation which 
is a typical coupled system, the methodology developed in this paper is potentially 
applicable to other coupled systems and more complicated port-Hamiltonian system- 
s. Some other important control subjects like the convergence of the minimal energy 
controls that drive the system state to rest in finite time in exact controllability and 
other control problems discussed in [10] and [30] are worthy being investigated in the 
future. 
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